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Abstract 



Tight-binding molecular dynamic simulations have revealed that Sii2 is an 
icosahedron with all atoms on the surface of an approximately 5 A diameter 
sphere. This is the most spherical cage structure for silicon clusters in the 
2-13 atom size range. 

PACS numbers: 36.40.+d, 61.43.Bn, 61.46,+w, 68.35.Bs, 82.65.My 
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Typeset using REVTgX 



The chemical reactions of silicon clusters with various reagents have puzzled both ex- 
perimentalists and theoreticians [[l -10]. Experimentally, Si^r clusters containing N = 2-50 



atoms have exhibited strong oscillations in reactivities as a function of cluster size |1|,@]. The- 
oreticians have attempted to explain this phenomenon either through the determination of 
cluster structures or through educated speculations PHT0| . Despite these efforts a consistent 
view of the cluster structures and reactivities has not emerged. 

Recently, we started to investigate this problem by calculating the structures of small 
silicon clusters using the tight-binding molecular dynamics (TB-MD) method developed by 
Menon and Subbaswamy [|TT|-|r4]]. This all valence electron method is equivalent to the ex- 
tended Hiickel method well known in theoretical chemistry fT5| . This method has yielded 



structures in excellent agreement with ab initio electronic structure calculations for both 
carbon and silicon clusters ||TT| p~4| . |T6]| . The particular Hamiltonian we employ consists of 



Harrison's universal tight-binding parameters ||17|| ) supplemented with two additional pa- 
rameters that are fit to the bond length and vibrational frequency of Si2 ■ Full details of 



this method are described elsewhere fl4| . Using this method in combination with classical 
molecular dynamics |18[ and a slow annealing schedule we optimized the structures of silicon 
clusters up to iV = 13. The cluster structures thus derived for sizes up to W = 11 are in 



complete agreement with the ab initio calculations of Raghavachari and co-workers [ I5 , Z0" |. 



This is remarkable considering that the parameters in this tight-binding Hamiltonian are not 
fit to any of the cluster structures ||14|| . This test verifies the reliability and predictability of 
our TB-MD method for the study of the structures of larger silicon clusters. 

We now focus our attention on N = 11-13 clusters because experiments with ethylene 
show the most dramatic variations in reactivities in this size range, with Sii2 being most 
reactive and Sii3 being least reactive, even though all these clusters are equally abundant 
0. We found that Sin is a tetragonal antiprism with one cap at the bottom and two caps at 
the top. Rohlfing and Raghavachari found this structure to be one of two possible candidates 
for the ground state structure of Sin [M. The structure of Sii3 may be described either as 
a 1-5-6-1 layered structure or as a capped trigonal antiprism, similar to that obtained by 
Rothlisberger and co-workers ||21|| , but unlike the atom centered icosahedron proposed by 
Chelikowsky ||22|| . Even though these two structures are interesting in themselves, neither of 



them are unique in their geometry. 

The TB-MD simulated annealing calculations also revealed that the lowest energy struc- 
ture of Sii2 is an icosahedron with a diameter of ~ 5 A. Distortion from the ideal icosahedral 
geometry is small. Such a highly spherical cage structure for this cluster had not been antic- 
ipated before. The cohesive energy of this structure is -3.75 eV/atom, its band gap is 1.8 eV, 
and its average coordination number is 5.0. Small distortions away from the ideal structure 
do not significantly change the cohesive energy. The ground state is a singlet in all cases, 
but the perfect icosahedron has a smaller band gap (0.26 eV) compared to its distorted form 
(1.8 eV). We have also investigated several other geometries for this cluster and found that 
they are about 0.4-0.7 eV higher in energy than the icosahedron. Hence, several structural 
isomers of Sii 2 are possible, with the icosahedron being the ground state structure. The local 
coordination of atoms in the icosahedral form of Sii2 is similar to that of the axial atoms in 
the pentagonal bipyramidal ground state of Sij. 

Our calculations also show that the cohesive energies of Sii2 and Sii3 are nearly the 
same, indicative that reactivity differences between these two clusters are not related to 
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their structural stabilities. However, Sin and Sii3 are either prolate or asymmetric tops, 
whereas Sii2 is a spherical top. This finding is consistent with the observation of Jarrold 
and co-workers that spherical clusters are at least an order of magnitude more reactive than 
the non-spherical clusters [p3[| . This is probably because the spherical surface provides the 
maximum number of nearest neighbor pairs of silicon atoms that are needed to bind the 
ethylene molecule. 

Since Sii2 is highly reactive it has not aroused as much theoretical interest as the inert 
clusters. Consequently, until now the structure of Sii 2 has not been investigated through 
accurate calculations. The few calculations P^J23| that do exist on Sii 2 obtained structures 
consisting of four- and five-membered rings that are very different from ours. However, these 
calculations employed the Stillinger- Weber |2j| or other empirical potentials that almost 
always yielded incorrect cluster structures, because these potentials were derived by fitting 
to the properties of bulk Si. Our tight-binding method does not suffer from this limitation 
because all the valence electrons are explicitly included in the calculation of the electronic 
energy 
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The icosahedral geometry with an atom at the center gives the closest packing of thir- 
teen hard spheres. The atoms interacting through two-body pair potentials typically form 
clusters with this structural pattern. For example, Ar 13 is an atom centered icosahedron 
|2Tfl . However, Ar 12 is not an icosahedron with an empty cage [57]. The metal clusters also 
prefer compact structures and do not allow empty cages |2l] , |28| . Thus the icosahedral cage 
structure is unique only for Sii 2 among elemental clusters. The compound clusters boranes 
and carboranes are the only other clusters that form icosahedral cage structures. 

If two-body interaction is dominant in our tight-binding Hamiltonian then we would ex- 
pect to obtain an atom centered icosahedron for Sii3. However, our calculations show such 
a structure to be about 2 eV higher in energy compared to the lowest energy structure. 
Furthermore, the structures of other Sijv clusters in the 3-11 size range are also not similar 
to those of Ar^ clusters. These two observations indicate that our tight-binding Hamilto- 
nian includes the necessary many-body interactions required to describe the silicon cluster 
structures correctly. 

In summary, our tight-binding molecular dynamic simulations predict icosahedral cage 
structure for Sii 2 . Such a highly symmetric non-crystalline cage structure has not been 
discovered previously for any other twelve atom elemental cluster. 

This research is supported by the New York University and the Donors of The Petroleum 
Research Fund (ACS-PRF $ 26488-G), administered by the American Chemical Society. 
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FIGURES 

FIG. 1. The lowest energy structure of the Sii2 cluster. It is a bi-capped pentagonal antiprism 
(BPA), called an icosahedron. This geometry is also describable as a 1-5-5-1 layered structure. 
Experiments have shown this cluster to be highly reactive [2]. 
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